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FOREWORD 


The present report is one of a series of three reports 
describing a new computer program which predicts turbulent 
boundary- layer behavior for a condition involving both heat 
transfer and pressure gradient. Part I serves as a summary 
report and describes the general analysis which is utilized in 
the numerical calculations scheme. In Part I I, the requisite 
low speed formulation consisting of a constant property flow 
with combined pressure gradient and mass transfer is described. 
Part III describes the numerical and computational procedures 
involved and serves as a computer program manual . 

The titles in the series are: 

PART I - Summary Report - "Calculation of Turbulent Boundary 
Layers with Heat Transfer and Pressure Gradient 
Utilizing a Compressibility Transformation," by 
C. Economos and J. Boccio. 

PART II - "Constant Property Turbulent Boundary-Layer Flow 
with Simultaneous Mass Transfer and Pressure 
Gradient," by j. Boccio and C. Economos. 

PART III - "Computer Program Manual," by J. Schneider and 
J. Boccio. 

This investigation was conducted for the Langley Research 
Center, National Aeronautics and Space Administration, under 
Contract No. NAS1-9624, with Mr. Kazimierz Czarnecki as the 
NASA Technical Monitor . 

The contractor's report number is GASL TR-748. 
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SUMMARY 


An analysis of the incompressible turbulent boundary layer, 
developing under the combined effects of mass transfer and 
pressure gradient, is presented in this paper. A strip-integral 
method is employed whereby two of the three governing equations 
are obtained by integrating the combined momentum and continuity 
equation to 50% and 100%, respectively, of the boundary -layer 
height. The latter equation is the usual momentum-integral 
equation; the former equation requires specification of shear 
at the point v=ri* = 0.5. Accordingly, Clauser ' s equilibrium 
eddy-viscosity law is assumed valid at this point. The third 
and final equation is obtained by specifying that Stevenson's 
velocity profiles apply throughout the domain of interest , from 
which a skin-friction law can be derived. 

Comparisons of the numerical results with the experiments 
of McQuaid, which include combined effects of variable pressure 
gradient and mass transfer, show good agreement. 
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CALCULATION. OF TURBULENT BOUNDARY LAYERS WITH 
HEAT TRANSFER AND PRESSURE GRADIENT UTILIZING A 
COMPRESSIBILITY TRANSFORMATION 


PART II - CONSTANT PROPERTY TURBULENT BOUNDARY - 
LAYER FLOW WITH SIMULTANEOUS MASS TRANSFER AND 
PRESSURE GRADIENT 

By J. Boccio and C. Economos 
General Applied Science Laboratories, Inc. 


INTRODUCTION 


In recent years there has been an increasingly greater 
interest in turbulent boundary layers. On one hand, continuing 
research has been aimed toward finding an understanding of tur- 
bulent shear flows of which the boundary layer is probably the 
most interesting example. On the other hand, an insistent demand 
still exists for reliable prediction methods to calculate stream- 
wise development of gross properties of the boundary layer as 
well as its mean velocity and shear-stress distribution. 

The serious interest in turbulent boundary layers has been 
evidenced by the fact that relatively large symposiums were 
dedicated to the subject. In 1968, two meetings were held - both 
of which served largely pragmatic interests. One held at Stanford 
University (Reference 1) dealt with the study of incompressible 
turbulent boundary layers, while the study of its compressible 
counterpart was the primary concern at the NASA sponsored meeting 
held at Langley Field, Virginia (Reference 2) . The 1968 Stanford 
Conference revealed many interesting aspects of the technical 
problem associated with the study of incompressible turbulent 
flows over impermeable surfaces subjected to both zero and non- 
zero pressure gradients. One outcome was the fact that some 
integral methods, specifically the dissipation-integral method 
of Alber (Reference 3) and the strip-integral method of Moses 
(Reference 4) were just as successful in predicting the mandatory 
test cases posed as the more elaborate differential or mean-field 
methods of Mellor and Herring (Reference 5), Cebeci and Smith 
(Reference 6) and Beckwith and Bushnell (Reference 7). 


At the Langley conference# the published theoretical pre- 
dictions by Economos (Reference 8) dealing with the study of 
compressible turbulent boundary layers over flat plates with 
mass addition and, more recently, the theoretical predictions 
made by Economos and Boccio (Reference 9) that treat the problem 
of compressible boundary-layer development with pressure gradient 
and heat transfer but no mass-transfer effects by utilizing a 
Coles' -type (Reference 10) compressibility transformation has 
revived, to some extent, the use of such types of transformations 
in the study of turbulent boundary-layer flows. 

As discussed in Part I of this report, reexamination of this 
compressibility transformation in an effort to improve velocity 
profile representation, has led to the development of a somewhat 
more general form of transformation. This new form relates the 
compressible boundary-layer flow under the influence of pressure 
gradient and heat transfer to an incompressible one in which 
pressure gradient and mass transfer occur simultaneously. Thus, 
in order to achieve the objectives of this contract, the first 
step was development of a formulation which described the latter 
flow. This report will describe such a formulation. 

An integral method of approach has been chosen because of 
its virtue of incorporating these effects on the turbulence 
structure in an implicit and global manner. Thus, the avoidance 
of local turbulence assumptions that are required in field 
methods, the fact that there is a dearth of corroborative experi- 
mental information which relates the dependence of mass transfer 
on eddy viscosity, and the simplicity of the method, all go into 
making the integral approach an attractive calculation tool for 
this endeavor. 

Accordingly, at the outset, a two-parameter profile repre- 
sentation is deemed necessary. And because of its firmer theo- 
retical basis, the Stevenson defect law (Reference 11) has been 
chosen over those postulated by Mickley and Smith (Reference 12) 
and McQuaid (Reference 13). The analysis could have used the 
defect law posed by Simpson (Reference 14) which is considered 
to be a logical extension to Coles "law of the wall - law of the 
wake formulation" to account for the non-constancy of the 
von Karman parameter, k^ , which occurs at lower momentum thick- 
ness Reynolds number (< 6000). but its added complication and only 
slight improvement did not mandate its use. From the defect 
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velocity representation emerges a skin-friction law from which 
a differential equation can be obtained relating the requisite 
three dependent variables that define the system; namely, the 
Coles wake parameter, n , the wall skin-friction parameter, o , 
and the Reynolds number based upon boundary-layer height, R . 
Consequently, in addition to this equation and the usual momen- 
tum-integral equation, a third differential equation is required 
to close the system. In this context, several approaches are 
available based upon their success reported at the 1968 Stanford 
conference. For example, higher moment equations such as the 
moment of momentum equation as Alber uses for the non-transpired 
problem can be utilized. And indeed, this has been the case in 
the method of Lubard and Fernandez (Reference 15) who have 
treated the identical problem as that presented herein, However 
to implement this method they must extend the procedure of Alber 
by not only uncoupling the attendant dissipation integral from 
the local pressure gradient, but they must also decouple its 
dependence from the mass transfer rate. Consequently, in addi- 
tion to using Alber' s empirical curve fit which relates the 
Clauser shape parameter with pressure gradient parameter and 
which they have to generalize to account for mass transfer, 
they must further assume special forms for the variation of mass 
transfer rate to associate with an "equilibrium" flow in order 
to close their system. Although the numerical examples cited 
are found to be weakly dependent upon which variation (three in 
all) is used and although they present rather good agreement 
with experiments, it is felt that use of such an assumption 
could be restrictive insofar as incorporating it into other 
types of pressure gradient and mass transfer rate distributions. 

Accordingly, the approach taken herein to obtain the remain 
ing governing equation and which is felt to be one degree less 
empirical with regard to its extension to problems involving 
transpiration is the strip-integral method of Moses - a method 
which has also enjoyed reasonable success at the Stanford con- 
ference. This additional equation is obtained by integrating 
the momentum equation up to 50% of the boundary-layer height, 
at which point the shear stress is a priori assumed to be that 
postulated by Clauser (Reference 16) for equilibrium flows. 

Thus, an assumption based upon equilibrium flows extended to 
account for mass transfer and arbitrary pressure gradient is 
required here as well as in the method of Reference 15. However 
this assumption is only required at one point in the boundary 
layer and not across the entire boundary layer as is required 



in the aforementioned reference. Also, no empirical curve fit 
of non-transpired bound ary- layer data a priori extended to the 
transpired problem need be made. 

In the following sections, the integral method is developed 
for the two-dimensional case and for the case where the injectant 
is the same fluid as the external stream. Numerical results are 
obtained for three specific problems, namely, those termed by 
McQuaid (Reference 17) as pressure distribution I, II, III. 
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SYMBOLS 


A. 

i 


elements of ordinary differential equations, 
c.f. Appendix and Eq. (38) 


A (.(f), v^) functional form, c.f. Eq. (24) 

B(v ) ( 1/4 ) v 

w w 

c. column vector relating to right-hand-side of the system 

^ of ordinary differential equations* c.f. Appendix and 

Equation (38) 


C f 

F i ( ^7*) 

G(ff,7j) 


local skin friction coefficient, 2 t /p u 2 

w e e 

i = 1,2,3 functional forms, c.f., Eqs . ( A-8 ) , ( A-9 ) , ( A- 10 ) 
functional form, c.f. Eq. (23) 


H 


form factor 


1 ( 77 , p, q) definite integrals described by Eq. (A-4) 
J(ir) functional form, c.f. Eq. (26) 


K 

k l 

k 2 

Pe 

Q i (77 , ,T7) 

R 

R- 

y 


- u/v 

constant, 0.4115 

constant, 4.9 

external pressure, lb/ft 3 

i = 1,...,5 functional forms, c.f. Eq. (A-5) 

Reynolds number based upon boundary-layer height, 

(p u ‘ 6/p) 

e, o 

Reynolds_number based upon local value of y 

(P u y/p) 
e , o 


U generalized velocity function, c.f. Eq. (8) 
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i >1 l«0 lo 10) i fa) Ip- Ip 


u 


u /u 
e e,o 


__ - ~ _ _ j, 

u ,u shear velocities, (r/.o) 2 , (r /p ) ~ , respectively, ft/sec 

T rw w 


u, v 

-+ -+ 
u , v 


V 

w 

w(tj) 


velocities respectively along and normal to body, ft/sec 

velocity ratios u/u ; v /u , respectively 
J Tw w rw J 

v /u 
w e,o 

difference between outer edge value and local value 
of Coles wake function, c.f. Eq. (19) 


w (rj) 


Coles wake function 


w. ( 17 ) ; i = 1,2,3 various curve fits to Coles wake function, 
1 c.f. Eq. (21) 


x,y 


-+ 

y 


s 
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space coordinates respectively along and normal to body 

pS Tw J/m 

value of y at laminar sub-layer height 

experimental constant (0.018) in the description of 
eddy viscosity c.f. Eq. ( 6 ) 

boundary layer height, ft 

displacement thickness, ft 

eddy viscosity, ft 2 /sec 

viscosity ratio, pe/u 

y/5 

a particular value of 7) taken within to be equal to 0.5 
momentum thickness, ft 

definite integrals, c.f. Eq. (31), (30), respectively 
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I 


jU molecular viscosity, lb sec/ft* 

v kinematic viscosity, ft 2 /sec 

ir Coles wake parameter 

p density, lb sec^/ft 4 

S 6/6 

r shear stress, lb/ft 2 

<p skin-friction parameter, u e /u^ w 

X Reynolds number defined as u (x-x ) /v 

e,o o 

6*/6 

Subscripts 

( ) pertains to local condition at edge of boundary layer 

( ) pertains to initial conditions 

( ) pertains to laminar sub-layer height 

s 

( ) pertains to condition at wall 

w 

Superscripts 

( ) normalization with respect to_ external value unless 

otherwise noted, i.e., u = u/u 

e 

( )* designates properties at 77* 

(t?*) 

( ) ' designates integration performed between 

the couple (0,rj*) 

( ) designates integration performed between 
couple (0,1) 

( ) 1 d/dx 
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ANALYSIS 


A. Basic Equations 

The differential equations describing the mean flow 
of an incompressible turbulent boundary layer for a two- 
dimensional flow field are taken to be 


Continuity: 


du dv 
Sx + 9y 


0 


Momentum : 


5u 


du 

W ay - 


d p e it 

dx + 


Shear Stress: 


r 


— . 3u 

^ 97 


u’V' = (jS+P«)|= 


External Flow: 


dp du 

e _ _ e 

dx ~°e e d 3T 


( 1 ) 

(2) 

( 3 ) 


( 4 ) 


The developing flow-field is subject to the boundary conditions 

u = 0 v = v (x) at y = 0 
w 

u = u (x) at y -* oo 

e 

A set of ordinary differential equations can be obtained 
fran the boundary- layer equations (1) and (2) by integrating 
in y yielding after some manipulations 
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-4 J^U’ d5 - S* £ J’’*udri - V -f- * 3dn) 


dX 


dX 


d\ 


+ ^ Un U e ){J‘ T? 2u a dtj - u* J V udrj - rj*) 


1 f~* 

I {t 




( 5 ) 


where the superscript ( ) denotes, unless otherwise noted, that 
the quantity has been ratioed with the local edge velocity, u .* 

In addition, the following definitions are required, namely: e 

* 5 G e,o ( *-*o )/f7 

n s y/6 


R = u 6/v 
e,o 

U = u /u 

e e e,o 

<p = u/(f/p ) i5 

e w 

f H f/p U= 

e e 

( ) s initial quantities 

( ) * s properties at i) a j)* 


( ) s external conditions 
e 

( ) s wall conditions 
w 


* To be consistent with Parts I and III of this contractual 
requirement, superscript ( ) connotes properties in the 
incompressible plane. 


9 


Any number of differential equations can be_obtained by eval- 
uating Equation (5) at different values of 77 *. Evaluation at 
rj* * 1 results in the usual momentum -integral equation. - 
Accordingly, with two unknown parameters assumed for the velocity 
profile representation (which is to be discussed later) , only one 
additional equation and choice of fj* is required. Values of 77 * 
between 0.2 to 0.5 have been tried by Moses and by Baron ti 
(Reference 18) with little effect on the final results. Thus, the 
arguments posed by these two authors are extended here by consider- 
ing as the second governing equation that which results by letting 
rj* = 0.5, at which point Clauser's eddy viscosity law is also 
assumed to hold, i.e. 


e* = 0 u 6* 
e 


0.018 u 6* 
e 


( 6 ) 


Considering the molecular viscosity negligible at 77 * = 0.5 
then the shear stress becomes 

r* - & f|>* - » O (#* < 7 > 

where the normalized displacement thickness is 

O = 6 */5 

B. Velocity Profile Representation 

The turbulent boundary layer is generally assumed to consist 
of three distinct flow regions: ( 1 ) the viscous sublayer, ( 2 ) the 
turbulent core, and (3) the wake region. A fourth region termed 
as a buffer or blending region which is bracketed by regions ( 1 ) 
and ( 2 ) above must also be considered in the description of a 
turbulent boundary layer. However, buffer region expressions are 
not as well agreed on and normally any expression for this region 
incorporates the. sublayer or the sublayer and turbulent core 
(Kleinstein, Reference 19). 


10 



In the viscous sublayer molecular viscosity dominates. in 
the buffer region, transition to the fully turbulent part of the 
boundary layer takes place and the laminar and turbulent mechanism 
appear to be of equal importance. In the turbulent core, molecular 
viscosity is assumed to completely lose its significance and the 
velocity profile depends only on the turbulence. Beyond the 
turbulent core is the wake region which is essentially another 
transition region in that the turbulence becomes intermittent and 
the flow changes from fully turbulent to the free stream condition. 

In incompressible flow over impermeable flat plates, various 
"laws" have been used to describe these regions. Regions (1), (2) 
and the buffer region, which is not considered any further in this 
report, are referred to as the "law of the wall" region. The wake 
region, so called because of the flow similarity to that of a 
wake propagating into an otherwise undisturbed flow, has been 
described by the "law of the wake" devised by Coles (Reference 20). 

The extension of these laws to describe the velocity profiles 
over permeable bodies are discussed in this section. The formula- 
tior is brief and represents a summary of existing work due to 
several contributors? additional discussions on this subject can 
be found in Reference 21. 


A three-layer model of the boundary layer is considered and 
within this framework, the choice of the profiles is believed to 
be the best representation available. 


A generalized velocity function defined as 

u = r u+ du*/u* 

" T 

is introduced where 

u + = u/u = ucp 

TW 

U - (r /P )* 5 

TW W 


u + = u /u 

r T TW 


( 8 ) 
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Now, in considering the ramifications associated with the 
definitions of the laminar sublayer and the law of the wall, 
namely, that within these regions the streamwise gradients do 
not significantly affect the shear distribution, reguires then 
that this generalized velocity function be 

U = (2/v + ) [ (l+vV’)* 2 - 1] (9) 

w w 


Furthermore, within the laminar sublayer where it is assumed 
that the quantity e , i.e., (pe/p) is considered much less than 
unity, it follows that the velocity distribution u of the form 


-+-+ 
u v 

w 


exp 



1 


0 < 


1 * 

y 


< 


- + 



do) 

where a measure of the laminar sublayer height, y , is, as yet, 
unknown and where 

— i- — ~ - 

v = v /u = v <p 
w w rw w 


-+ 

y 


3 pu r u = 

TW 


R U rj/<p 


With regard to the law-of-the-wall region where e 
considered to be much larger than unity then 


-+ 


is 


-+ 

e 




( 11 ) 


Accordingly, by considering the momentum transport theory of 
Prandtl with its attendant mixing-length hypothesis, namely that 

e + = k.u + y + (12) 

X T 

yields, when substituted into the above equation and integrated, 
U = k 2 + (l/k 1 ) In y + y + S y^ (13) 
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The constants k_ and k^ are taken to be those associated with 
the zero mass-transfercase. A posteriori justification for 
this assumption has been provided by experimental measurements 
for the zero pressure-gradient case with mass injection and is 
discussed in Reference 21. For pressure-gradient cases its 
justification will be provided by comparisons with the 
experiments. 

Finally, by combining Equation (9) with Equation (13) there is 
obtained the law-of-the-wall equation with transpiration in 
the form 


U=[(l+v^u + ) is -l] = k 7 + (lA,) In y + ,• y^ y* 

w “ ’ (14) 

which is identical to that deduced by Stevenson. To complete 
tljie profile formulation within these two regions, the value of 
y must be determined. Requiring the velocity to be continuous 
at this point requires the equality of Equation (10) and 
Equation (14) which is provided if 

l+(v*/2) [kj+tlA^n y^] = exp [ ( 1/2 ) v^ y * ] (15) 


With the velocity profiles determined within the inner 
turbulent region, the variation in the outer region, the 
so-called defect law, is readily derivable fraxi a dimensional 
similarity argument. Thus, by considering the existence of a 
region of common validity between the law-of-the-wall and the 
law-of-the-wake and the fact that the latter law avoids a 
direct confrontation with the physical mechanism of shear 
turbulence requires the defect law to have the form 

f(u* , v* ) - f(u + , v+ ) = F 1 (rf) (16) 

Identifying the function, f, with the generalized 
velocity function, U, then requires that F ( 77 ) must be 

F (rj) = ( 2 /v* ) [(L+uV )*-(l+uV (17) 

W G W W 
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and, by considering the universality of F(tj), then its repre- 
sentation can b$ taken from that which is associated with the 
impermeable wall problems. Thus, in the spirit of Coles, F (f>) 
behaves like 

F(r>) = (ir/k x ) W(n) - (lA^ in rj ( 18 ) 


where 


W (77) = w(l) - w (77) 


( 19 ) 


and wf?7 ) is the so-called Coles Wake function. Thus, the defect 
law is obtained by equating Equation ( 17 ) with Equation ( 18 ) , 
thereby yielding the general representation within the outer 
layer , i . e . > 

U = ( 2 /vt )[(l+vt u*)* 5 -!] = k 2 + (lAJ-tny* + 

W W Z L 

(n A-^wfn) 

which further demonstrates the utility of the generalized 
velocity function U. 

Thus, the velocity profile representation, taken in 
concert, is 

Inner Region, law-of- the -wall region 

u + v + = exp (v + y + ) -1 0 < y ^ y* (20a) 

w • w s 

(2^ >C(l«V; )**-!] = y + ?*> F* (20b) 

* 

where y must satisfy the trancendental equation 
s 

1 + (v + /2) [k r) +(lA 1 ) -t-ny* ]=exp[.(l/2)v + y* ] (20c) 

w z 1 s w s 
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I 


Outer Region, law-of-the-wake region 


( 2 /v + )[(l+v + u + ) ^** 1 3 = (lAJ^n y + + (ff/k. )w(t 7 ) + *2 

W W 1 ± 


P < V * 1 
' s 

y + <y + *J + 

2 s 2 

Various curve fits to the Coles wake function have been 
incorporated into the numerical program, namely 

w (rj) = w 1 (t 7 ) = 2 


w (v) = w 2 ( 77 ) = l-cosirn (21) 

w(n) = w 3 ( 77 ) = 3pTy 3 -125f7 4 +l83iT 5 -13 3rj 6 +38r) 7 

However, only the results using w^ ( 77 ) are presented herein. 

A working form of the velocity -pro file representation is 
obtained from Equation (20d) by normalizing the velocity 
component with TT e a nd the physical height with 6 yielding 
the following velocity- profile representation, which is con- 
sidered to be more suitable for integration and differentiation. 

Thus, the velocity profile is given by 

u = 1+A ( 0 , v^) G (ff, rf) +B (v^) G (ir,rj) (22) 

where 


G (ff • 7?) 

® [J(17)W(t?)+K 

1 — 1 

fP 

c 

(23) 

A(^,v w ) 

2 h 

s -(Vv 

(1 AS) 

(24) 

B <v w ). • 

d/4) V w 


(25) 

J (ff) = 

(irAij) 


(26) 

K H - (1A X ) 


(27) 


Note, Equation (22) is strictly applicable outside the laminar 
sublayer region; thus this region is neglected when this equation 
is substituted in the integrals associated with Equation (5) . 
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C. The Skin Friction Law 

A direct consequence of Equation (20d) is the skin-friction 
law obtained by placing rj equal to unity at 7 = 7. Slight 
manipulations yield 

(1+^5) =1+ (v*/2) [kj+flA^) -t-n 6 + +(irA 1 )w(l) ] (28a) 


which can be rewritten as 

(l+^ 2 v )** = l+(v w ^p/ 2 ) [k^flA-^ ^n(RU e /^) + ( 2 ffA 1 ) ] ( 28 b) 

with, according to the curve fits of the wake function, 
w(l) =2 

This equation, which relates the skin-friction parameter, 
a scaled boundary- layer thickness, IT, and the wake parameter, 
it, to the independent variable, x, constitutes the remaining 
equation of the problem at hand. Implicit in the equation is 
the x-wise dependency of v^ since and u are considered to 
be known functions of x. For later computational convenience 
differentiating Equation (28b) with respect to x yields after 
considerable manipulations a working differential equation 
namely 

(-tnR) ' (lA 1 )+(-^nU e ) ' (lA 1 ) + (ff) ' (2A l )+(^) ' {-(1A X ) + 

+ (2/^ w )[(lVv w ) !j ’ 1] ^(lVv/V^Vw) '{ (2/®v w ) C(1 + 

p 2 ? ) Ji -l]-^(l-^ 2 v )"*} = 0 (29) 


where 

( ) ' = d/dx 

It can easily be shown that in the limit of v -*0 the above 
equation reduces to that for the impermeable case which 
comprises part of the analysis found in Reference 9. 
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D. Final Working Forms of the Governing Equations 

By examining the form of the velocity profile given by 
Equation (22) and the integration required in Equation (5), it 
is deemed more convenient in defining integrals of the form 


nT^ 

3* 

ill 

T? 

/ 

O 

(u-1) 2 d tj 

= A2 (^,^ w ,ir,i7) 

(30) 

(7?) 

A i a 

f 

o 

(u-l)d^ = 

Ai (<p,v w ,ir,r)) 

(31) 


where it can be readily shown that the term multiplying 
d/dx(^nR), in the latter equation can be replaced with 

+ (2-u*) A 1 ( ^* ) + (1-u *)rf* 

Now, by considering the functional dependency of these integral 
forms on the dependent variables, 75> ff, R and their implicit 
dependence on the independent variable, j£, through the possible 
streamwise variation of v and u , then the requisite 
differentiation of these Yntegra? forms when substituted 
back into Equation (5) will, after collecting terms, produce 
the result that 

(trip) + (it)' F^ 77 * ) + (-tnR) ' [A^ t? * ) + (2 — u*)a| 7? * ) + 

(l-Q*) n* 1 + (-t-nU ) ' { [2A^ r? * ) + (4-u*)A,^* ) + (1-u*)??*] - 
_ e _ it G* 1 

F< v * } } + (InVj 'F 2 (77 * } + ( -=—)-( = )It*~ 2 2 ] = 0 

9 

(32) 


where V = v /u 

w w e,o 

Noting that the auxiliary equation described above will 
yield identically the so-called momentum-integral equation if 
the integration is performed up to 77 = 1 (with u*and r* equal 
to one and zero respectively at this point), then it can be shown, 
formally at least, that the last and final governing equation 
can be written as 
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I 


(Irq) 'F + (ir) -F 3 (1) + (-triE) ' [A^ 1 ^ + A X U) ] 


( 1 ) 


+ (^nu e ) ' [2 aJ 1) +3A^[ 1) - F^ 1} ] + (^nV w ) *F 2 
+ (tf /^) + (l/R(p 2 ) = 0 


( 1 ) 


(33) 


where the superscript (1) indicates that the functions are to 
be evaluated at f? - 1* Also, in terms of the definitions, 

Ajl and Aj, the normalized momentum and displacement thickness 
are respectively 


L = 0/6 = - [ A2 (1> + A 1 <1> 3 


O = I */6 = - A x (1) 


(34) 


(35) 
(r?*> (**> 


The .necessary equations tg evaluate the terms F * , . 2 

F <*>. and the terms A -^ 17 , A 2 ^ » are given in the 

Appendix^ Also, in lieu of the above equations, the shear 
stress, t*, can now be formally represented by 


r* " mbG/brj) [A(p,tf )+2 b(v )G(ff,r»)] (36) 

w w 

where 


(3G/5i7)=[j (ff) +K/?f] 


(37) 


Finally, Equation (29) , (32) ,and (33) can be recast into 
the form 


A 1 A 4 A 7 


( trup ) ' 


c i 

A 2 A 5 A 8 


(-tdR) 1 

= 

(N 

U 

A 3 A 6 A 9 


( n ) ' 


C 3 


where the coefficients A^ and C are also given in the Appendix. 
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The integration of the system (38) can be readily carried 
out by means of a high-speed computer, with the development of 
the boundary layer thug being determined once the external 
pressure field, the mass-transfer distribution, and the initial 
conditions are prescribed. 

E. Initial and Boundary Conditions 

For boundary conditions all that is required is a specifi- 
cation of the external and wall velocity distribution and the 
unit Reynolds number. The initial conditions which are 
necessary in order to numerically integrate the above system 
of equations are the values of Jp, it, and R at some initial station, 

V 

Usually, integration will start at a station where there i 
an experimentally determined velocity profile. Then a plot of 
versus Ry can be generated from which R can be ascertained. Also, 
by making several judicious guesses on jp, a 'log-law', i.e.. 
Equation (14) which best describes the data can be ascribed. 

The result is a generalized 'Clauser Plot', as shown in Figure 1 
from which the desired value of Jp can be obtained. With ip and 
R now known, the third parameter, i r, is obtained by solving 
the skin-friction law. Equation (28b) > which is recast 
below in a slightly different form, i.e., 

ff= (k x /2) { [ (l-Hp 2 ^) h -l] (v^Jp/2 ) ~ 1 --k 2 -(l/k 1 ) In (RU & fip) } (39) 

where by definition U q is unity. 

If one wished to initiate the solution at a leading edge, 
the assumption implicit in the analysis are that at x q =0 

C^) Q = 13.55 

(it) o = 0.6 

from which the skin- friction law, yields an expression for 
R , namely 

R = exp[ (0 ,06073/v )[(1 +183.603v ) 2 -l]-0.6l04} (40) 

w w 

for the leading edge value once a (^) is prescribed. 
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WIXI 


Data, Ref (17) , Pressure Dist 
= 8.3 in 

= 7247, v =.00826 

w 



Eq. (14) 


nnc 




RESULTS AND DISCUSSIONS 


The numerical results of the analysis of the previous 
section have been compared with the experimental data of 
McQuaid (Reference 17) . In the three cases considered, the 
pressure gradient and blowing rate were varied simultaneously 
but were restricted to air into air injection without heat 
transfer. Boundary- layer developments were measured with in- 
jection rates (v /u ) between 0 and 0.0008 at freestream 
velocities of 50 W an§'°150 ft/sec. At the outset, it must be 
stated that the skin-friction comparison is not made with the 
experimental values quoted by McQuaid. It is felt, and as 
McQuaid concedes, that considerable errors can accrue if use 
is made of the momentum-integral equation to experimentally infer 
the wall skin-friction value. Instead, -and as McQuaid suggests, 
the skin-friction coefficient used in all the examples cited 
is that which has been estimated by means of the Clauser Plot 
method previously described in Section E . In general, the 
agreement in all three cases is quite excellent. 

McQuaid Pressure Gradient I. - in this case, the boundary 
layer has developed under a mild adverse pressure gradient and 
a mild increase in blowing rate. Figure 2 shows the excellent 
agreement between the analysis and the skin-friction and 
integral data over the entire range. The ability of the 
analysis to predict the velocity profiles is also demonstrated 
in the next figure. in Figure 3, only three profile comparisons 
are shown; one at the forward portion of the region of interest, 
one in the center, and one at the end of the measurement domain. 
Ten such profiles have been measured by McQuaid and all have 
been compared with the existing analysis showing similar and 
ofttimes better agreement than those presented herein. 

McQuaid Pressure Distribution II . - A mild favorable 
pressure distribution, together with a decrease in the blowing 
rate, gives rise to the experimental data reported by McQuaid for 
this test configuration. Again, the agreement of skin friction, 
form factor, momentum thickness, and' velocity profiles as 
portrayed in Figures 4 and 5, is quite good. 

M cQuaid Pressure Distribution III - In this final case, the 
severest pressure gradient is considered. Over a span of approx- 
imately 30 inches the external velocity has increased by a 
factor of two while the blowing rate has increased three-fold. 
Again, the comparisons as exemplified in Figures 6 and 7 show 
rather remarkable agreement. Also, as before, all ten profiles 
reported have been compared, all showing similar comparative 
agreement . 
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FIGURE 5: TYPICAL VELOCITY PROFILE CCMPARISON,MCQUAlD DATA, PRESSURE DISTRIBUTION II 








— Theory 

A Experiment - McQuaid (17) 



FIGURE 6. COMPARISON OF TNOORV AND EXPERIMENT, 
PRESSURE DISTRIBUTION NO. til 





CONCLUSIONS 


An integral method has been developed for calculating the 
development of an incompressible boundary layer under the 
simultaneous influence of both mass transfer and pressure 
gradient. No need to generate new empirical constants has been 
required; and also there has been no need to utilize curve fits 
of non- transpired boundary layer data in the numerical program. 
The numerical program is quite simple and fast, taking approxi- 
mately 20 seconds to complete a case during which time ten 
theoretical velocity profiles can also be generated. No sub- 
stantial effort has been made to compare this analysis with the 
integral method of Lubard and Fernandez; suffice to say that 
both methods show similar comparative tendencies. However, it 
is felt that the approach present herein is less empirical than 
that of Lubard and Fernandez, or stated somewhat differently, 
the empiricism required manifests itself only at a discrete 
point of the boundary layer, namely, the assumption that 
Clauser’s eddy-viscosity law is valid at only one point, namely 
at fj * = 0.5. Hence, there is no need to assume gross profile 
behavior as is required in Reference 15, where 

1. a curve fit of Clauser's equilibrium form factor 
with pressure gradient is required. 

2. this same curve fit must be extended to include 
the transpired problems , and 

3 . a need to define an equilibrium flow with mass 
transfer is mandatory. 
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APPENDIX A 


I. The Profile parameters , A2 

Considering the velocity profile represented by Equation 
(22) and, for convenience rewritten here, i.e.. 


where 


u = 1 + A(<p,v w )G(ir,rj)+B(v w )G 


(A-l) 


G ( 7 T, 77 ) 5 [ J (IT) W (77) +Kln 77 J 
A (<p,v w ) = - (1+yp. v w ) 2 (l/(p) 


B (V ) - ( 1 / 4 ) v 

w w 


J(ff) = ( 1 /^) (IT) 

K = - (l/k 1 ) 

W (77) = w(l) - w ( 77) 


U 5 u/u , V = V /u 
e w we 


and substituting into the definitions for A^,A^ yields formally 
two expressions of the form 


A, (7?) = A(«,v )Q (ff,rj)+B(v )Q (IT, 77) (A— 2 ) 

1 ~ W 1 W Z 

A 2 ^ = A 2 (<p, V^) Q 3 (ff, 77) + 2 A (^,V w )B(V w )Q 4 (TT,77)+B 2 (vMC^ ( 7T , f? ) 

(A-3 ) 


where the superscript, (n) implies the 77-dependency of the 
functions brought about by the integration interval (0,77) . 
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The expressions for the Q's involve integrals in the form 


I (t), p,q) 


7? 

= J* [W (fj) ] P (ln77) q d 77 . 

o 


(A-4 ) 


where the powers p and q represent integers ranging from 0 to 
4. Accordingly it can be shown that 


Qi ( 7 T, T 7 ) = J (it) I(rj, 1,0) + KI ( 77 , 0 , 1 ) 

Q 2 (TT,T7) = J 2 (ff) I (77,2,0)+2j(ff)KI (ij, 1,1)+K 2 I(T7,0,2) 

Q 3 (ff,f 7 ) = Q 2 (tT,77) (a-5) 

Q 4 (7T , 77 ) = J 3 (ff)I (77 / 3 / 0)+3J 2 (ff)KI(77,2,l)+3j(TT)K 2 I(^ / l,2) 

+k 3 i (77,0,3) 

q 5 (it ,ff) =j 4 (tt) i ( 77 ”, 4,0) +4j 3 (tt)ki(t?„ 3 ,o)+ 6 j 2 (it) k 2 i ( 77 , 2 , 2 ) 
+4J (tt) K 3 ! (77, 1. 3 ) +K 4 I (fT, 0 , 4 ) 


and specification of a curve fit to the Coles wake function, 
i.e. , 

w ( 77 ) = 2 (3 77 2 - 2^) (A- 6 a) 


or 


w (ff) = 1 - cos (TTfj) (A-6b) 

or 

W(f?)=39 n 3 - 125^ 4 + 183fj 5 - 133r) 6 + 38T7 7 (A-6c) 

is all that, is required to evaluate the Q's for any r). 
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II. Derivatives of the Profile Parameters 

Examining the auxiliary equation in the body of the report, 
i.e.. Equation ( 5 ) indicates that these profile parameters 
must also be differentiated since the term 


V* 2 * 7?* 

(d/dX) J u drf - u(d/dx) j udjj 
o ° 


can be written as 


V* 2 *?* 

(d/dx) / [(u-1) + 2u-l]dr/ - u* (d/dx) J C(u-l)+l]dr7 
o o 


and since 77 * is a constant (taken here to be 0.5) then the 
original term in the auxiliary equation can be replaced by 

(d/dx) A 2 ( ^* ) + (2-3)*(d/dX) A 1 (f? * ) 


Now, considering the functional dependency of the A’ s on 
!o, It, v and that 

~ W 



(X) 




u 

e, o 



o 


u 

e 


y* > 

V*’ 


Then the above expression when differentiated and common 
terms collected can be expressed symbolically as 


(d/dx) A 2 (f ’* ) + (2-u)*(d/dx)A 1 ( ^* ) = 

(ln^p) 1 (F l ) (f ’* ) + (lnV w ) ' (F 2 ) (7? * ) + (ff) 1 (F 3 ) (r) * ] 
- (lnU e ) 1 (F 2 ) (71 * ) 


(A-7) 
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where the expressions for the F's are 


F n (fr * ) = (dA/d <p) (7p) [2AQ 3 +2BQ 4 +(2-u)*Q 1 3 (A-8) 

F, ( ^* } = (dA/dv ) (v ) [2AQ_+2BQ. + (2-u)*Q ] + 

Z W W 3 4 1 

(dB/dv J (v ) [2AQ .+2BQ, +(2-u)*Q ] (A-9) 

W W 4 D 1 

F^ A 2 (dQ 3 /dff) +2AB (dQ 4 /dtr) +B 2 (dQ 5 /dff) + 

(2-u)*[A(dQ 1 /dfr)-!-B(dQ 2 /dff) ] (A-10 ) 


Note, that for convenience, the indication of dependency of 
the various functions on the governing variables has been 
dropped here. Furthermore, it can readily be shown that 


(dQj/dir) = (dj/dff) K^rl»0)-(lA 1 )i(r?,l,0) 

(dQ 2 /dff) = (dJ/dir) [2Ji(rj,2,0)+2KI(77,l,l) ] 

(dQ 3 /dff) = (dQ 2 /d ir) (A- 11 ) 

(d0 4 /dff) = (dJ/dir) [3J 2 I (^,3,0) +6JKI (f?, 2 , 1 ) +3k 2 I (f?, 1 , 2 ) .] 
(dQ 5 /dff) = (dJ/dtr) [4j 3 l(77,4,0)+12j 2 KI(f?,3,l)+12JK 2 I(n,2,2) 
+ 4K 3 I (ff, 1,3) ] 
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Likewise the various derivatives of A are 


(dA/dp) = (l/<p 2 ) (l+^ 2 v ) 

w 

(dA/dv w ) = -(1/2) ty) (iVv^ 

and for B all that is required is 

(dB/dv ) = (1/4) 
w 
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Ill Elements of System (38) 

a ={-(lA.) + (2/^v ) Ul+^v ) is -l]-2 <p (i+^r 2 v )~ H ] 

-L *L W W W 


a 4 = {lA^ 

A 7 = {2A 1 ) 

C 2 = -{(inu e ) 1 (lA^) + (^nv w ) '{ (2/^v w ) [a^vj ^-1] 


Note: Since v = v /u = (v /u ) (u A ) = (V /U ) 
w w e w e,o e,o e we 

A = F ( ^*> 

2 1 


= <2-u>*A<”*> 


•5 *2 

A = F ( ^* ) 
A 8 3 


(2-u) a:" '+(i-a ) V 


= -{(-t.nU e ) , {[2A 2 ( ^* ) + (4-u)*A ; [ T? * ) + (l-u)*rj*]- F 2 ( ^* ) } 

__ V u* 

+ (tnV .) '{F 0 (T7 * } }+ if* - l/<p 2 }} 




A_ = F. 


( 1 ) 




A = F. 


( 1 ) 


C 3 = - {(^nUj , {2A 2 (1) +3A 1 (1) -F 2 (1) } + (toV w ) '{F 2 (1) } 


+ (V w /R) + (l/R)^p 2 j 
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